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We study the finite-size scaling (FSS) property of the correlation ratio, the ratio of the correlation 
functions with different distances. It is shown that the correlation ratio is a good estimator to 
determine the critical point of the second-order transition using the FSS analysis. The correlation 
ratio is especially useful for the analysis of the Kosterlitz-Thouless (KT) transition. We also present a 
generalized scheme of the probability-changing cluster algorithm, which has been recently developed 
by the present authors, based on the FSS property of the correlation ratio. We investigate the two- 
dimensional quantum XY model of spin 1/2 with this generalized scheme, obtaining the precise 
estimate of the KT transition temperature with less numerical effort. 
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The development of efficient Monte Carlo algorithms is 
important for studying many-body problems in physics. 
We recently developed a cluster algorithm, which is called 
the probability-changing cluster (PCC) algorithm, of lo- 
cating the critical point automatically [1]. It is an ex- 
tension of the Swendsen-Wang (SW) algorithm [2], but 
we change the probability of cluster update (essentially, 
the temperature) during the Monte Carlo process. Since 
we do not have to make simulations for several param- 
eters, we can extract information on critical phenomena 
with much less numerical effort. We applied the PCC 
algorithm to the 2D diluted Ising model [3], investigat- 
ing the crossover and self-averaging properties. We also 
extended the PCC algorithm to the problem of the vec- 
tor order parameter [4] with the use of Wolff's embedded 
cluster formalism [5]; studying the 2D classical XY and 
clock models, we showed that the PCC algorithm is also 
useful for the Kosterlitz-Thouless (KT) transition [6]. 

In the original formulation of the PCC algorithm [1], 
we used the cluster representation of the Ising model 
(generally, the Potts model) due to Kasteleyn and For- 
tuin [7] in two ways. First, we make a cluster flip as in 
the SW algorithm [2]. Second, we change the probabil- 
ity of connecting spins of the same type, p, depending 
on the observation whether clusters are percolating or 
not. We use the finite-size scaling (FSS) relation for the 
probability that the system percolates, E p , 

E v (3>,L) = X{tL}l v ), t=(p c -p)/ Pc , (1) 

to determine the critical point. Here, L is the system 
size, p c is the critical value of p for the infinite system, 
and v is the correlation-length critical exponent. We may 
alternatively consider that t = (T — T c )/T c , where T is 
the temperature. With a negative feedback mechanism, 
we locate the size-dependent temperature, E p of which is 
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1/2. The point is that E p has the FSS property with a 
single scaling variable. We may use quantities other than 
E p which have a similar FSS relation. Then, we could 
generalize the PCC algorithm for a problem where the 
mapping to the cluster formalism does not exist. 

In the FSS analysis of the simulation, we often use 
the Binder ratio [8] , which is essentially the ratio of the 
moments of the order parameter m. The moment ratio 
has the FSS property with a single scaling variable, 

T% = f^l (2) 

as far as the corrections to FSS are negligible. Here, the 
angular brackets indicate the thermal average. However, 
the moment ratio derived from a snapshot spin config- 
uration is always one; therefore, the instantaneous mo- 
ment ratio cannot be used for the criterion of judgment 
whether we increase or decrease the temperature. 

In this Letter, we treat the correlation ratio, the ra- 
tio of the correlation functions with different distances. 
Studying the thermodynamic properties of the correla- 
tion ratio, we will show that it is a good estimator for 
determining the critical points of the second-order and 
KT transitions. We also propose a generalized scheme of 
the PCC algorithm based on the correlation ratio. Using 
this general scheme combined with the quantum Monte 
Carlo simulation of the continuous-time loop algorithm, 
we study the 2D quantum XY model of spin 1/2. 

Let us start with considering the spin-spin correlation 
for the D dimensional systems, 

9(r) = ^J2 Si - Si + r ' ( 3 ) 

i 

where Si is a spin at site i, and N = L D is the number of 
spins. Precisely, the distance r is a vector, but we have 
used a simplified notation. We have assumed a transla- 
tional invariance. For an infinite system at the critical 
point, the correlation function decays as a power of r, 

<g{r))~r-<- D -**i\ (L = oo, t = 0), (4) 
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with the decay exponent r\. Away from the critical point, 
the ratio of the distance r and the correlation length £ 
plays a role in the scaling of the correlation function. 
Moreover, for finite systems, two length ratios come in 
the scaling form of the correlation function; 

(g(r, t, L)) ~ r-^- 2 ^ h(r/L, L/£), {L + ^t + 0). 

(5) 

Then, the ratio of the correlation functions with different 
distances r and r' becomes 

(g(r,t,L)) _ ( r_y{D-2+r,) h(r/L,L/£) 
,r'J 



(9(r',t,L)) 



h(r>/L,L/ty 



(6) 



If we fix two ratios, r/L and r/r', the correlation ratio 
takes the FSS form with a single scaling variable, 



(9(r,t,L)) 
(9(r',t,L)) 



(7) 



In case the correlation length diverges with a power law, 
£ oc t~ v , Eq. (7) becomes the same form as Eq. (2); how- 
ever, we should note that Eq. (7) is also applicable to the 
case of the KT transition, where the correlation length 
diverges more strongly than the power-law divergence. 

In order to examine the FSS properties of the correla- 
tion ratio, we simulate the 2D Ising model on the square 
lattice with the periodic boundary conditions. As for the 
distances r and r' , we choose L/2 and L/4; we take the 
horizontal or vertical direction of the lattice for the orien- 
tation of two sites. We show the temperature dependence 
of the moment ratio (m 4 ) / {m 2 } 2 and the correlation ratio 
{g{L/2))/{g{L/A)) in Fig. 1. From now on, we represent 
the temperature in units of J/ks, where J is the cou- 
pling constant and ks is the Boltzmann constant. The 
error bars are the order of the width of curves. We see 
from Fig. 1 that the crossing of the data of different sizes 
are better for the correlation ratio; in other words, the 
corrections to FSS for the correlation ratio are smaller 
than those for the moment ratio. We check the two ra- 
tios at the critical point for a quantitative comparison. 
The size dependence of the moment ratio for the 2D Ising 
model was carefully studied by Salas and Sokal [9]. The 
correlation functions of the critical 2D Ising model for 
L — > oo can be calculated with the use of the continuum 
field theory [10]. The moment ratio is written in terms 
of the integral of the 2n-point correlation functions, and 
the fourth-order moment ratio becomes 



(m 4 )/(m 2 y = 1.1679229(47) 



(8) 



in the limit L — > oo [9] , where the numbers in the paren- 
theses denote the uncertainty in the last digits. Without 
performing an integral, the correlation ratio is simply 
written as 

(g(L/2)) |<M1/2)|-V4 £^^(1/4)1 



< 5 (L/4)> 1^(1/4)1-1/4 £4 =1 |0„(l/8)|' 



(9) 



where 9 v {z) are the Jacobi ^-functions. Thus, we can 
easily obtain the correlation ratio at the critical point as 
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FIG. 1: Plot of the moment ratio (m 4 )/(m 2 ) 2 and the corre- 
lation ratio {g(L/2))/{g(L/4)) of the 2D Ising model for L = 
8, 16, 32, and 64. The critical values of two ratios for L = oo 
are given by open circles. 



in the limit L — > oo. The critical values of two ratios for 
L = oo, Eqs. (8) and (10), are shown by open circles in 
Fig. 1. The calculated correlation ratio at T c = 2/ ln(l + 
V2) = 2.269- • • for L = 8 is 0.9418(5). The deviation 
from the value for the infinite lattice is only 0.2% even 
for L = 8, which is smaller than that for the moment 
ratio, 0.7%. We mention that the scaling plot of the 
correlation ratio as a function of (T — T C )L 1 ^ V with v = 1 
is very good. Our results suggest that the correlation 
ratio is a good estimator to determine T c using the FSS 
analysis. 

Let us turn to the KT transition. As an example, we 
pick up the 2D 6-state clock model. It was shown that 
the 2D g-state clock model has two phase transitions of 
the KT type at T x and T 2 (Ti < T 2 ) for q > 4 [11]. We 
here simulate the 2D 6-state clock model on the square 
lattice. We show the temperature dependence of the mo- 
ment and correlation ratios in Fig. 2. For the correlation 
ratio, the curves of different sizes overlap in the inter- 
mediate KT phase (T\ < T < T 2 ), and spray out for 
the low-temperature ordered and high-temperature dis- 
ordered phases, which is expected from the FSS form of 
Eq. (7). Then, we can make a FSS analysis based on the 
KT form of the correlation length, 



£ oc exp(c/\/i). 



(11) 



(g(L/2))/{g(L/4)) = 0.943904982 



(10) 



We try scaling plots of (g(L/2))/(g(L/A)} as a function 
ofL/exp(ci/Vl -T/Ti) and L/ cxp(c 2 / yfT/T 2 - 1) for 
the low- and high-temperature sides, respectively. Using 
the data for L=8, 12, 16, 24, 32, 48, and 64, we esti- 
mate two KT transition temperatures. The best-fitted 
estimates arc Ti=0.6980(22) and T 2 =0.9008(14), which 
are compatible with the recent results using the PCC al- 
gorithm [4], Ti=0.7014(ll) and T 2 =0.9008(6). On the 
contrary, as seen from Fig. 2, the corrections to FSS are 
larger for the moment ratio, which makes the FSS analy- 
sis difficult. We have shown that the correlation ratio is a 
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FIG. 2: Plot of the moment ratio (m 4 )/(m 2 ) 2 and the corre- 
lation ratio (g(L/2))/{g(L/4)} of the 2D 6-state clock model 
for L = 8, 16, 32, and 64. 



FIG. 3: Plot of Tkt(L) of the 2D quantum XY model with 
5=1/2 for L = 8, 12, 16, 24, 32, 48, 64, 96, and 128, where 
I = InfeL. 



good estimator especially for the KT transition. It is due 
to the fact that we only use the property of correlation 
function; the characteristic of the KT transition is that 
the correlation function shows a power-law decay at all 
the temperatures of the KT phase. We will discuss this 
point later. 

We may use the FSS properties of the correlation ra- 
tio for the generalization of the PCC algorithm. In- 
stead of checking whether the clusters are percolating or 
not, we ask whether the instantaneous correlation ratio 
g(L/2)/ '<?(L/4) is larger or smaller than some fixed value 
R c . Of course, we can use other sets of distances. We 
decrease (increase) the temperature, if g(L /2) / 'g(L / '4) is 
smaller (larger) than R c . We start the simulation with 
some temperature. We make the amount of the change of 
temperature, AT, smaller during the simulation; in the 
limit of AT — > 0, the system approaches the canonical 
ensemble. In most cases near the critical temperature, 
g(L/2) and <?(L/4) take some positive values for ferro- 
magnetic systems. It sometimes happens that g(L/2) 
and/or g{L/A) becomes or negative. We simply regard 
such a case that the system is disordered, and then in- 
crease the temperature. 

Here, we apply the generalized scheme of the PCC al- 
gorithm to the study of the quantum spin system. The 
loop algorithm, especially the continuous-time loop algo- 
rithm [12], has been successfully used for various quan- 
tum spin systems. We deal with the 2D S = 1/2 quan- 
tum XY model as an example, and determine the critical 
point automatically. The Hamiltonian is written as 



H = -J^iSfSJ + S?S?), 



(12) 



where the spin operators S x ' y are one half of the Pauli 
matrices a x - y . We are concerned with the correlation 



(13) 



which is an off-diagonal element if we take z axis as a 
quantized axis. Although off-diagonal elements arc diffi- 
cult to calculate with the conventional quantum Monte 
Carlo simulation methods, they are calculated easily with 
the continuous-time loop algorithm [13]. We should note 
that when calculating the correlation with the improved 
estimator, the correlation becomes non-negative. 

We treat the systems with linear sizes L = 8, 12, 16, 
24, 32, 48, 64, 96, and 128. We start with AT = 0.005, 
and gradually decrease AT to the final value, 0.0001. Af- 
ter 20,000 Monte Carlo sweeps of determining T C (L), we 
make 10,000 Monte Carlo sweeps to take thermal aver- 
age; we make 100 runs for each size to get better statis- 
tics and to evaluate the statistical errors. We calculate 
g{L/2)/g(L/A) to check whether it is larger than R c or 
not; the value of R c is set to be 0.8. 

We follow the same procedure for analyzing the KT 
transition as the 2D classical XY model [4]. Using the 
FSS form of the correlation ratio and KT form of the cor- 
relation length (Eqs. (7) and (11)), we have the relation 



Tkt(T) — Tkt + 



(InfeL) 



2 ' 



(14) 



We plot Tkt(T) as a function of l~ 2 with I = InfeL for 
the best-fitted parameters in Fig. 3. The error bars are 
smaller than the size of marks. Our estimate of Tkt is 
0.340(1). We estimate the uncertainty by the x 2 test of 
the data for 100 samples. This value is compatible with 
the estimate of the recent study, 0.3427(2) [14]. The 
constant c, in Eq. (14), is estimated as c = 2.45(5). 

To discuss the critical exponent rj, let us consider 
the correlation (g(L/2)) at Tkt(L). In Fig. 4, we plot 
(<?(L/2)) as a function of L in a logarithmic scale. In 
order to estimate 77, we use the FSS form including small 
multiplicative logarithmic corrections 



(,g(L/2)) = AL- r '(lnb'Ly 



(15) 



the existence of the multiplicative logarithmic corrections 
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FIG. 4: Logarithmic plot of (g(L/2)) at T K t{L) of the 2D 
quantum XY model with 5=1/2 for L = 8, 12, 16, 24, 32, 48, 
64, 96, and 128. 

was pointed out for the KT transition [15, 16]. We should 
note that the logarithmic corrections are negligible for 
the correlation ratio. Using the above form, we obtain r\ 
= 0.251(5) and r = 0.055(9). We show the fitting curve 
obtained by using Eq. (15) in Fig. 4. This value of r\ is 
consistent with the theoretical prediction, 1/4 (=0.25). 
Our logarithmic-correction exponent r takes almost the 
same value as that for the classical XY model; 0.038(5) 
[4] and 0.0560(17) [16]. 

It should be noted that Alet and S0rensen [17] also 
applied the PCC algorithm to the quantum Monte Carlo 
simulation. One checks the ratio of the correlation length 
and the system size in their method. 

In summarizing, we have shown that the correlation 
ratio is a good estimator to determine the critical points 
of both the second-order and the KT transitions. A gen- 
eralized scheme of the PCC algorithm based on the FSS 
property of the correlation ratio has been also presented. 

Although the study of the correlation ratio was partly 
motivated by a possible application to an efficient Monte 



Carlo algorithm, the FSS analysis of the correlation ratio 
itself is of great interest. The moment ratio [8] has been 
employed in the study of various problems of the phase 
transition, from the lattice gauge theory to the polymer 
statistics. The use of the correlation ratio may help pre- 
cise studies of the critical properties for complicated sys- 
tems. The FSS analysis for the KT transition sometimes 
faces a difficulty due to the singular divergence of the 
correlation length and the logarithmic corrections [18], 
and may lead to false conclusions. As a way to escape 
from such a difficulty, the level-spectroscopy method [19] 
was developed for treating quantum systems which show 
the KT transition. Nomura [18] studied the rcnormaliza- 
tion of the correlation functions for the 2D sine-Gordon 
model, and elucidated the logarithmic corrections on the 
fixed lines. It can be shown that the leading corrections 
are canceled if we consider the ratio of the correlation. 
This is the reason why the correlation ratio is a very 
good estimator especially for the KT transition, which 
has not been noticed before. In contrast, the corrections 
may not be canceled for the moment ratio. In Fig. 2 
we have given a direct demonstration of the intermediate 
KT phase of the 2D clock model with rather small sizes, 
which clearly indicates the superiority of the correlation 
ratio in the analysis of the KT transition. The FSS of 
the correlation ratio is very promising for the study of 
various systems which show the KT transition. 

Our generalized scheme of the PCC algorithm is now 
free from the restriction of the cluster mapping. It can 
be applied to many problems. For example, the cluster 
formalism does not work well for frustrated systems, but 
we can use the generalized scheme of the PCC algorithm. 
We can also apply the generalized scheme to a wide va- 
riety of quantum systems. 
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